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ABSTRACT 

We calculate Lyman Alpha Emitter (LAE) angular correlation functions (ACFs) at 
z — 6.6 and the fraction of lifetime (for the 100 Myrs preceding z ~ 6.6) galaxies spend 
as Lyman Break Galaxies (LBGs) or as LBGs with Lyman Alpha (Lya) emission using 
a model that combines SPH cosmological simulations (GADGET-2), dust attenuation 
and a radiative transfer code (pCRASH). The ACFs are a powerful tool that sig¬ 
nificantly narrows the 3D parameter space allowed by LAE Lya and UV luminosity 
functions (LFs) alone. With this work, we simultaneously constrain the escape fraction 
of ionizing photons f esc = 0.05 — 0.5, the mean fraction of neutral hydrogen in the 
intergalactic medium (IGM) (xhi) ~ 0.01 and the dust-dependent ratio of the escape 
fractions of Lya and UV continuum photons f a /fc = 0.6 — 1.2. Our results show that 
reionization has the largest impact on the amplitude of the ACFs, and its imprints 
are clearly distinguishable from those of f esc and / Q // c . We also show that galaxies 
with a critical stellar mass of M* = 10 8 5 (10 9 5 )Mq produce enough luminosity to 
stay visible as LBGs (LAEs). Finally, the fraction of time during the past 100 Myrs 
prior to z = 6.6 a galaxy spends as a LBG or as a LBG with Lya emission increases 
with the UV magnitude (and the stellar mass M*): considering observed (dust and 
IGM attenuated) luminosities, the fraction of time a galaxy spends as a LBG (LAE) 
increases from 65% to 100% (~ 0 — 100%) as Muv decreases from Mjjv = —18.0 to 
—23.5 (M* increases from 10® — 10 1o ' 5 Mq). Thus in our model the brightest (most 
massive) LBGs most often show Lya emission. 

Key words: radiative transfer - methods: numerical - dust, extinction - galaxies: 
high-redshift - dark ages, reionization, first stars 


1 INTRODUCTION 

The Epoch of Reionization (EoR) marks a major phase 
change in the ionization state of the Universe. While the 
intergalactic medium (IGM) is predominantly composed of 
neutral Hydrogen (H i ) at the beginning of this epoch, it 
is completely ionized by the end, as a result of H i ionizing 
photons produced by both stars and quasars. However, the 
progress of reionization has been hard to pin down since 
it depends on a number of parameters including the initial 
mass function (IMF) of reionization sources, their star for¬ 
mation rates (SFR), their stellar metallicity and age, the 
escape fraction of H i ionizing photons produced by each 
source and the clumping of the intergalactic medium (IGM), 
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to name a few. A further complication is introduced by su¬ 
pernova feedback and (to a lesser extent) the ultraviolet 
background (UVB) built up during reionization in suppress¬ 
ing the gas content (and hence star formation) in low-mass 
galaxies which are the main sources of H i ionizing photons 
(see e.g. Barkana fe Loeb|2001 Ciardi & 


et al.|20lT Sobacchi fe Mesinger|2013||Wyithe fe Loeb|2013| 


IMaio 


Dayal et al.||2015 and references therein). 

Lyman Alpha (Lya) photons are a powerful tool in un¬ 
derstanding the ionization state of the IGM given their high 
optical depth (r) to H i (e.g. Madau fc Rees|2000) 


t = i.5 x io 5 h 


1 + z 


3/2 


(1 + S H ) Xhi, (1) 


where h is the Hubble parameter, Uli, and represent the 
cosmic baryon and matter density, respectively, (1 + 5h) is 
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the hydrogen over-density and \hi is the fraction of neutral 
hydrogen. As seen from this equation, even a neutral hydro¬ 
gen fraction as low as 10 _o can lead to a significant atten¬ 
uation of Lya photons at higli- 2 , making them extremely 
sensitive probes of H i in the IGM. As a result, a class of 
high -2 galaxies called Lyman Alpha Emitters (LAEs), de¬ 
tected by means of their Lya line (at 1216 A in the rest 
frame), have become popular probes of reionization, with 
statistically significant samples available in the reionization 


epoch, at 2 ~ 5.7 and 6.6 ( 

Malhotra et al.| 

2005 

Taniguchi 

et al. 

2005 

Shimasaku et al. 

2006 Hu et al. 2010 

Oashikawa 

et al. 

2006 

Ouchi et al. 2010 

Kashikawa et al. 201 

l). Indeed, 


a number of theoretical papers have used semi-analytic (e.g. 


Dijkstra et al. 2007 Dayal et al. 2008 Samui et al.||2009 

or numerical (e.g. 

McQuinn et al. 2007 

Iliev et al. 201IX 

Dayal et al. 2011 

Forero-Romero et al.| 

2011, Duval et al. 

2014 Hutter et al. 

2014) models to reproduce the observed 


number counts of LAEs (the Lya luminosity function; LF) 
at various redshifts in the reionization epoch. However in¬ 
terpreting a change in the Lya LF is rendered challenging 
by the fact that the observed Lya luminosity depends on: 
(a) the fraction of H i ionizing photons (l-/ eS c ) produced by 
a galaxy that are able to ionize the interstellar H i resulting 
in the Lya recombination line, with the rest (/ eBC ) escaping 
to ionize the IGM, (b) the fraction of the intrinsic Lya pho¬ 
tons that can emerge out of the galactic environment unat¬ 
tenuated by dust (/□,) and H i and (c) the fraction of these 
emergent Lya photons that are transmitted (T Q ) through 
the IGM (depending on ( Xhi )) and reach the observer. 


Given these uncertainties, we require an alternative 
measurement to constrain the ionization state using LAEs. 
One such strong measure is provided by the two-point an¬ 
gular correlation function (ACF) of LAEs that describes 
their spatial clustering. Indeed, [McQuinn et al.| (2007) have 
shown that the spatial clustering can hardly be attributed 
to anything other than the large scale (~ Mpc scale) ion¬ 
ization regions created during reionization. This result is 
similar to that obtained by Jensen et al. (2014) who find 
that upcoming large-field LAE surveys should be able to 
detect the clustering boost for sufficiently high global IGM 


neutral fractions (20% at 2 = 6.5), although Jensen et al. 


120131 point out that a sample of several thousand objects 
is needed to obtain a significant clustering signal. Further¬ 
more, Zheng et al. (2010) have shown that Lya radiative 


transfer modifies the ratio of observed and intrinsic Lya lu¬ 
minosities depending on the density and velocity structure 
of the environment, i.e. LOS (transverse) density fluctua¬ 
tions are suppressed (enhanced), leading to a change in the 
amplitude of the two-point correlation function compared 
to the case without the environmental selection effect. How¬ 
ever, while Behrens & Niemeyer (2013) have confirmed a 
correlation between the observed Lya luminosity and the 
underlying density and velocity field, they do not find a sig¬ 
nificant deformation of the two-point correlation function by 
post-processing hydrodynamical simulations with a Lya ra¬ 
diative transfer code. 


A second probe is presented by the fraction of Ly¬ 
man Break Galaxies (LBGs) that show Lya emission: given 
that the physical properties of LBGs do not evolve in the 


150 Myrs between 2 ~ 6 and 7, a sudden change in the 
fraction of LBGs showing Lya emission could be attributed 


to reionization (e.g. Stark et al.||2010j [2011 Fontana et al. 

2010 

Pentericci et al. 2011 

Caruana et al. 2014; Faisst et al. 

2014 

Schenker et al. 

2014 

Tilvi et al. 2014 

|. However, this 


interpretation is complicated by caveats including the red- 
shift dependence of the relative effects of dust on Lya and 
ultraviolet (UV) continuum photons (Dayal fe Ferrara|2012 1 


and the fact that simple cuts in EW and UV luminosity may 


lead to uncertainties in the LAE number densities (Dijkstra 
& W yithe|2 012 I. 


Coupling a cosmological SPH simulation snapshot at 
2 ~ 6.6 with a radiative transfer (RT) code (pCRASH) and 


a dust model, Hutter et al. (2014) have shown that the ef¬ 
fects of /esc , fa and T a are degenerate on the LAE visibil¬ 
ity: reproducing the observed Lya LFs can not differenti¬ 
ate between a Universe which is either completely ionized 
or half neutral (( Xhi }— 0.5 — 10 —4 ), or has an / eac ranging 
between 5-50%, or has dust that it either clumped or ho¬ 
mogeneously distributed in the interstellar medium (ISM) 
with f a /fc = 0.6 — 1.8; here f c represents the fraction of 
UV photons that emerge out of the ISM unattenuated by 
dust. In this work, we extend our calculations to use the 
two-point LAE ACF to narrow down the 3-dimensional pa¬ 
rameter space of (xhi) , /esc and fa/fc■ Further, we study 
the time evolution of both the Lya and UV luminosities to 
show the fraction of time during the past 100 Myrs prior to 
2 = 6.6 that a galaxy is visible as a LBG or as a LBG with 
Lya emission, both considering the intrinsic and observed 
(dust and IGM-attenuated) luminosities. In principle, the 
ratio of these fractions could be measured by relating the 
number of observed LBGs with Lya emission to the number 
of observed LBGs at z = 6.6 — 7.3. 

The cosmological model corresponds to the ACDM Uni¬ 
verse with dark matter (DM), dark energy and baryonic den¬ 
sity parameter values of (Da, D m , Qb) = (0.73, 0.27, 0.047), 
a Hubble constant Hq = lOO/i = 70kms _1 Mpc _1 , and a 
normalisation as = 0.82, consistent with the results from 
WMAP5 (Komatsu et al. [2009 ). 


2 THE MODEL 


In this section we briefly describe our physical model for 
high-redshift LAEs that couples cosmological SPH simu¬ 
lations run using GADGET-2 with a RT code (pCRASH, 
Parti et al. (20111) and a dust model, and interested readers 
are referred to Hutter et al. (20141 for a detailed description. 


The hydrodynamical simulation analysed in this work 
was carried out with the TreePM-SPH code GADGET-2 
and has a box size of 80/U 1 comoving Mpc (cMpc) and 
contains 1024 3 DM particles, and initially the same num¬ 
ber of gas particles; the mass of a DM and gas particle 
is 3.6 x 1O 7 1 i _ 1 M0 and 6.3 x 10 6 /i _1 Mq, respectively. The 
simulation includes all the standard processes of star forma¬ 
tion and its associated metal production and feedback using 


the prescription of Springel & Hernquist (20031, assuming 
a Salpeter ( Salpeter||1955 ) initial mass function (IMF) be¬ 
tween 0.1 — lOOM®. Bound structures of more than 20 to- 
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^esc -0-05 f esc -0.25 f esc -0.50 



1. Angular correlation function for simulated LAEs. The mean ACF is calculated from 36 mock catalogues (12 along each of 
Erections) assuming volumes with a depth of 30/i —1 Mpc and a FoV of ~ 3 x 10 3 h —2 Mpc 2 . In each panel, the solid line shows the 
the best fit ( fa./fc , fesc , ( Xhi )) combinations for which the simulated LAE Lya LFs are within the lcr limit of the observations 
likawa et al.| ( |2011| with shaded regions showing the variance across the mock catalogues. In each panel, the grey dashed line 
le ACF for LBGs from the whole box, corresponding errors on the ACF for LBGs are comparable with the line width, and black 
epresent the observational results by |Kashikawa et al.| ( (2QQ6| ). Columns show the results for f eS c = 0.05,0.25 and 0.5, as marked. 
Lies of (xhi) are marked at the end of each row, with the fa/fc value marked in each panel, along with the x 2 error. 
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tal (Dark Matter, gas and star) particles are recognised as 


galaxies using the Amiga Halo Finder (AHF; Knollmann & 


Knebe 2009). Of all these galaxies, we only use “resolved” 


galaxies that are complete in the halo mass function in all 
our calculations - these consist of at least 160 (10) gas 
(star) particles and have a halo mass Mh XL 10 9,2 Mq. As¬ 
suming each star particle to have formed in a burst, we 
calculate its spectra, and rest-frame intrinsic Lya (L(, ni ) 
and UV continuum (L 1 ™/.; 1505A) luminosities depending 
on the stellar mass, age and metallicity using the popula¬ 
tion synthesis code STARBURST99 ( jLeitherer et al.|1999 |. 
For each galaxy, the dust mass and its corresponding UV 
attenuation are calculated using the dust model described 
in Dayal et al. (2011). The observed UV specific luminosity 
is then calculated as L ob f. = f c x L 1 ™/. where f c is the frac¬ 
tion of UV photons that emerge out of the ISM unatten¬ 
uated by dust which is fixed by matching to the observed 
evolving UV LF at z ~ 6 — 8; galaxies with an absolute 
UV magnitude Muv ^ — 17 are then identified as LBGs. 
The observed Lya luminosity (L/f s ) is then calculated as 
L* = f a T a L where f a and T a account for dust attenu¬ 
ation in the ISM and by H i in the IGM, respectively. Galax¬ 
ies with L° a ba ^ 10 42 erg s" 1 and EW = L obs / Lff c ^ 20 A 
are identified as LAEs. 

Given that both L™ 1 and (xhi) depend on the frac¬ 
tion of H i ionizing photons that can escape out of the ISM 
(/e S c), we use five values of / eac =0.05, 0.25, 0.5,0.75, 0.95 
to post-process the z ~ 6.6 snapshot of the hydrodynam- 
ical simulation with the 3D radiative transfer (pCRASH). 
To explore the full range of (xhi) (that determines T a ) we 
run pCRASH starting from a completely neutral to a com¬ 
pletely ionized state and obtain ionization fields for different 
values of the mean fraction of neutral hydrogen (xhi) ; for 
each galaxy we assume the initial line profile to be Gaussian 
and compute the average Lya transmissions along 48 differ¬ 
ent lines of sight (LOS). For each combination of f esc and 
( xhi) the transmission T a is then fixed and the only free 
parameter that can be adjusted to match Lya luminosities 
to observations is the relative escape fraction of Lya and 
UV photons from the ISM, p = f a // c . 

From the Lya LFs we find the following trends: the 
amplitude of the LAE Lya LF decreases with increasing 
f eac due to a decrease in L/ff*. However, this decline can be 
compensated by an increase in either f a /f c or T a (in an 
increasingly ionized IGM). Comparing our model Lya LFs 
to observations (Kashikawa et al. 20111 we found that al¬ 
lowing for clumped dust (p A. 0.7), the observations can be 
reproduced for ( xhi)— 0.5 — 10 —4 , f esc — 0.05 — 0.5 and 
fa/fc = 0.6 — 1.8 within a Ict error. 


3 CONSTRAINTS FROM LAE CLUSTERING 

Given the sensitivity of Lya photons to H i , the ionization 
field should be imprinted in the spatial distribution of the 
LAEs which can be quantified by the two point correlation 
functions (ACF). In this section, we use the ACF of theoreti¬ 
cal LAEs for each of these best-fit combinations (that match 
the Lya LF to within a 1 -a error as summarised in Table [l]) 


Table 1 . For the ( XHI) value shown in Column 1, we summarise 
the f a /f c ratio required to best fit the simulated Lya LF to ob¬ 
servations (jKashikawa ct al. 2(11 ij within la limits in columns 2-4 
with subscripts showing the f esc value used (|Hutter et al.|20 1. 


(Xhi) 

(fa/fc) 0.05 

(fa/fc) 0.25 

{fa/fc) 0.5 

0.50 

1.2 

1.6 

1.8 

0.25 

0.8 

1.2 

1.4 

0.10 

0.68 

1.0 

1.4 

0.01 

0.68 

0.9 

1.2 

10~ 4 

0.60 

0.8 

1.2 


to narrow the joint constraints on ( Xhi ), fesc and fa/fc ■ We 
start by describing our procedure for obtaining LAE ACFs, 


which depends both on the depth along the LOS (Peebles 


1980) since all galaxies are projected to a plane perpendic¬ 


ular to the LOS, as well as the chosen field of view (FoV), 


which should be comparable to that observed (Kashikawa 
et al.|2006). Indeed, while the measured ACF should be in¬ 


dependent of the FoV if large enough areas are sampled, the 
restricted observational FoV leads to an ACF that is not 
independent of sample variance. To get an estimate of the 
average ACF and its variance, we generate 36 mock cata¬ 
logues (12 along each of x,y,z directions) for over-lapping 
volumes comparable to that observed by |Kashikawa et ah] 
12006), corresponding to a redshift distance Az ~ 0.1 at 
z ~ 6.5 and a field of view (FoV) of ~ 3 x 10 3 h _2 Mpc 2 . 
Using the Landy-Szalay estimator we compute the ACF in 
each mock survey region denoted by Wi(r) with respect to 
the mean LAE number density n of the complete simula¬ 
tion box and estimate the mean value of w(r) as well as its 
variance from our mock catalogues as 


(1 + w(r)) = ^5Z ni ( 1 + Wi ( r )) 


( 2 ) 


We start by calculating the LBG ACF (over the entire 
box) to get an estimate of the underlying galaxy popula¬ 
tion. As seen from Fig. |T| LBGs are almost homogeneously 
distributed and the ACF is consistent with essentially no 
clustering on scales ,$ 30 h~ Mpc. On the other hand, the 
LAE ACF is affected both by ISM dust, as well as the large 
scale topology of reionization. It might be expected that in 
the early stages of reionization, only those galaxies that are 
clustered and hence capable of building large H ii regions 
would be visible as LAEs (leading to a large amplitude of 
the ACF), with the amplitude of the ACF decreasing as 
reionization progresses and faint objects are able to trans¬ 
mit enough flux through the IGM to be visible as LAEs. 
Indeed, as shown in Fig. [l] LAEs exhibit precisely this be¬ 
haviour. For a given / esc value (we remind the reader this 
corresponds to a fixed L ((“*), as the IGM becomes more ion¬ 
ized (going down the vertical columns in the Fig. [Tj) , smaller 
galaxies are able to transmit more of their flux through the 
IGM, requiring lowering f a /fc values to fit the Lya LF. 
While for a half neutral IGM (panel a), only strongly clus¬ 
tered galaxies are visible as LAEs ( w(r ) ~ 4.5) at scales of 
~ 2/U 1 Mpc, T a increases for a completely ionized IGM at 
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Figure 2. The grey shaded regions show the (1 — 5 )o (dark grey to white respectively) regions for the combinations of fesc, ( Xhi ) and 
fa/fc that best fit the observed Lya LF data fronTHutter et al. ( 201 4 ). We ove r-plot black contours for 3 o (stars), 5 a (dots), > bo 
(hatching) by comparing model results to observed ACF data |Kashi kawa_et al. ( 2006 ). As shown, LAE clustering observations (3cr) 
require an IGM that has \HI ~ 0 . 01 , fa/fc '1 1.2 and fesc -C 0 . 5 . See Sec.|3|for details. 


the same scale, leading to a more homogeneous LAE distri¬ 
bution resulting in a lower amplitude of the ACF (panel m; 
w(r) ~ 1.5). 

At a given value of (xhi) , L „ decreases with increas¬ 
ing fesc (horizontal rows in the same figure) which must be 
compensated by an increase in fa/fc- However, this compen¬ 
sation results in very similar ACFs at a given ( \hi) value. 
Our results therefore show that the ACF is driven by the 
reionization topology (as determined by ( Xhi )), with / esc 
and the local f a /fc having a marginal effect (a factor of 
about 1.5) on its amplitude. 

We then calculate the yf errors between our simulated 
ACFs and observations, and find that observations con¬ 
strain ( xhi) ~ 0.01 for f esc = 0.05, 0.5 and ( xhi) A 10 -4 
for fesc = 0.25 (to within a 3cr error). Further, while the 
fa/fc ratio is compatible with homogeneously distributed 
dust for fesc = 0.05, the decrease in L™ 1 requires clumped 
dust ( fa/fc ^ 0.7) for fesc = 0.25 and 0.5. 

To highlight the importance of the spatial clustering of 
LAEs, we show the (1 — 5)cr constraints allowed by matching 
the Lya LFs to observations in Fig. [2] As seen, these encom¬ 
pass a region such that ( xhi)— 10“ — 0.5, fesc = 0.05 — 0.5 
and fa/fc = 0.6 — 1.8. However, building ACFs for each of 
these allowed combinations, we find that theory and obser¬ 
vations yield much tighter constraints of {xhi)— 0.01 —10~ 4 , 
fe a c = 0.05 — 0.5 and f a /f c = 0.6 — 1.2 to within a 3<r error. 

Finally, our results show that it is the reionization 
topology (as parameterised by {xhi)) that drives the ACF, 
supporting the results obtained by McQuinn et al. (2007), 


Jensen et al. (2013) and Jensen et al. (20141. Although 


Jensen et al. (2014) have assumed a simple scaling down 


of the L a luminosity of all galaxies by a fixed amount due 
to galaxy evolution, they also find that galaxies are more 
likely to be observed as LAEs if they reside in ionized re¬ 


gions for ( xhi) ~ 20%. We note that the number of galaxies 
we identify as LAEs is similar to within 10% for the different 
fesc, (x«r) an d fa/fc combinations, demonstrating that the 
enhanced LAE clustering in our model can be attributed to 
an increasing neutral IGM. 

We note that the average number of LAEs in our mock 
catalogues (~ 300) exceeds the number of identified ob¬ 
jects in Kashikawa et al. (20061. In order to compare to 
a complete sample we have considered all galaxies with 
L°f s ^ 10 42 erg s _1 and EW ^ 20A as LAEs. Since this 
luminosity cut may include fainter galaxies than observed, 
the obtained ACFs represent a lower limit, as the cluster¬ 
ing increases for higher luminosity cuts ( Jensen et al.|2014 1. 
According to Jensen et al. (20131 our galaxy sample size of 
~ 300 is sufficient to distinguish a half-ionized from an ion¬ 
ized IGM, but does not provide the necessary sample size of 
~ 500 to distinguish ( Xhi)= 0.25 from ( Xhi)= 1- However, 
we have modelled the luminosity of each galaxy according 
to its stellar population obtained from the hydrodynami- 
cal simulation. This makes our results more sensitive to the 
ionization state of the IGM than Jensen et al. (20131; their 
clustering signal was reduced by the random scatter they 
added to the imposed fixed mass-to-light ratio. Hence, in 
agreement with McQuinn et al. (2007) we find our sample 
size sufficient to distinguish {xhi)= 0.25 from {xhi)= 10 -4 . 


4 THE RELATION BETWEEN LAES AND 
LBGS 

In this section we show how simulated z ~ 6.6 galaxies 
build up their stellar mass (AL*), and the time evolution of 
their UV and Lya luminosities, and Lya equivalent widths 
(EW). We then present the fraction of time during the last 
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Figure 3. Time evolution of M*, Lya: and UV luminosities and EWs across three M* bins of 10 s , 10 9 and 1O 1U M0 as marked above 
each column. The upper 3 panels show the intrinsic quantities (M* in black; L™* in red (upper line); LV lt in blue (lower line)); we use 
/esc = 0 for the Lyo luminosity. The middle three panels show the observed quantities (M* in black; L^ s in red (upper line); in 

blue (lower line)) where we assume T a = 0.45, f a = 0.68 f c and a individual / c depending on the dust mass of each galaxy. The dashed 
lines in the top two panels show the current observational limits corresponding to the Lya (10 42 erg s -1 ) and UV (10 39 ’ 6 erg s -1 /A) 
luminosities. The lower 3 panels show the intrinsic (yellow, upper line) and observed (green, lower line) Lyo EWs with the dashed line 
showing the minimum limit of 20A. In each panel shaded regions show the variance in the given M* bin. 


100 Myrs prior to z = 6.6 for which galaxies of different UV 
magnitudes are visible as LBGs or as LBGs with detectable 
Lya emission, both considering intrinsic and observed lu¬ 
minosities for each combination of / esc , (xhi) and f a /fc 
that reproduces the Lya LF and LAE ACF within lcr and 
3<J respectively. In order to consider a time span indepen¬ 
dent of resolution effects we have chosen a period of time of 
100 Myrs as an interval. We discuss the limiting case, i.e. 
considering the total lifetime of each galaxy in our simula¬ 


tion as the corresponding period of time, in the Appendix 


4.1 Time evolution of stellar mass, Lya and UV 
luminosities 

We use the ages of each star particle to trace the growth of 
M*, and the metallicity and time dependent values of the 
intrinsic UV and Lya luminosities (L" l ‘ an d L ™ 1 , respec¬ 
tively) and the intrinsic Lya EW (= L™*/L l ™* c ) for galaxies 
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Figure 4. Fraction of time during the last 100 Myrs prior to z = 6.6 that galaxies spend as LBGs {fLBG: circles) and as LBGs with a 
Ly a line {f LBGa , squares), as a function of the intrinsic UV luminosity for intrinsic values of Lya and UV luminosities. The panels show 
the fractions for the indicated values of f eac = 0.05, 0.25, 0.5. The mean stellar mass in each Mjjy bin is encoded in the shown colour 
scale. The fractions are computed as the mean of the galaxies within Muv bins k ranging from k — 0.25 to k + 0.25 for k = —25 .. —18 
in steps of 0.5. Error bars show the standard deviations of the mean values. The dotted black line in each panel represents f LBGct f° r 
feac — 0.05; as clearly seen, increasing / esc to 0.5, leads to a decrease in /lbGq- 


in three bins of Af* ~ 1O 8,9 ’ 1 o M 0. We assume that all ion¬ 
izing photons emitted within galaxies are absorbed in the 
ISM and produce Lya radiation, i.e. f esc = 0, to investigate 
the most extreme case. 

As seen from the upper three panels of Fig.|3j both 
and L™ 4 (averaged over the galaxies in the given M* bin) 
rise as galaxies steadily build up in (stellar) mass, albeit 
with a large scatter reflecting the assembly history of differ¬ 
ent galaxies. As expected for normal star forming galaxies, 
the average intrinsic Lya EW for all the three M* bins con¬ 
sidered here has a value between 30 — 300A (lower 3 panels 
of the same figure) that is larger than the minimum value 
of 20A used to identify LAEs. Indeed, we find that once 
a galaxy exceeds a critical mass of roughly 10 8 ' 5 (10 7 ' 5 )Mq, 
it can produce enough luminosity to intrinsically be a LAE 
(LBG), it has also met the LBG criterion. 

As for observed luminosities, we remind the reader that 
we compute Lff s using the fraction of UV photons that es¬ 
cape out of the galaxy (/ c ); L°f s is computed assuming ho¬ 
mogeneously distributed dust (f a /fc = 0.68) and an IGM 
transmission value of T a = 0.45 (H utter et al.|2 014l. As ex¬ 
pected, including the effects of dust and the IGM reduces 
both the Lya and UV luminosities, as shown in the middle 
panels of the same figure so that the critical M* at which 
a galaxy has the minimum luminosity to be a LAE (LBG) 
increases to 1O 9 ' 5 (1O 8,5 )M0, although in most cases the ob¬ 
served EW is still larger than the minimum required value 
of 20A (bottom most panels). 

As expected from our discussion above, the time a 
galaxy spends as a LBG or LAE increases with increas¬ 
ing AL*. However, in our model a galaxy becomes a LBG 
before it also turns into a LAE because of the more strin¬ 


gent (luminosity + EW) constraints imposed on identifying 
a galaxy as a LAE, as expected from the critical M* val¬ 
ues quoted above. To quantify, while galaxies with stellar 
masses of ~ 10 s M@ are visible as LBGs for roughly the last 
80 Myrs, they do not meet the selection criterion to be vis¬ 
ible as LAEs: although their intrinsic EW are larger than 
the minimum value of 20 A required by observations (bot¬ 
tom panel of same figure), L)," 4 does not meet the required 
Lya luminosity of 10 42 erg s _1 . Galaxies with M* ~ 1O 9 M0 
are massive enough to sustain star formation and maintain 
and L™ 1 values above the required limits. Again, the 
lifetime spent as a LAE is driven by L“‘, with the intrinsic 
EW always exceeding 20 A(bottom panel of Fig. [ 3 ]). 

When considering observed luminosities, the time a 
galaxy is visible as a LBG decreases; the decrease is more 
pronounced for LAEs that are additionally affected by IGM 
transmission. From panels (d)-(f) of Fig. [ 3 } we see that in ad¬ 
dition to galaxies with 1 O 8 M 0 , galaxies with M* ~ 1 O 9 M 0 
are also no longer visible as LAEs as a result of the Lya lu¬ 
minosity dropping below visible limits. 


4.2 The fraction of time spent as LBG and as 
LBG with Lya emission 

We now calculate the fraction of time during the last 
100 Myrs prior to z = 6.6 that galaxies in different UV 
magnitude bins spend as a LBG ( fLBG ) and as a LBG with 
Lya emission (fLBGc i.e. as a LAE), and discuss the result¬ 
ing fractions (/ibgo and /ibg) for galaxies of varying UV 
magnitudes Muv = —18 to —25 in more detail. We start 
with intrinsic UV luminosities as shown in Fig. [4] firstly, 
as galaxies become more massive, they are able to sustain 
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Figure 5. Fraction of time during the last 100 Myrs prior to z = 
6.6 that galaxies spend as LBGs (/igg, circles) and as LBGs with 
a Ly a line (fLBGa, squares) as a function of the UV luminosity 
for our best fit models. The mean stellar mass in each Muv bin 
is encoded in the shown colour scale. This panel shows the best 
fit case for fesc = 0.05, ( Xhi)= 10 -4 and fa/fc = 0.60. The 
fractions are computed as the mean of the galaxies within Miry 
bins k ranging from k — 0.25 to i + 0.25 for k = —25 .. —18 
in steps of 0.5. Error bars show the standard deviations of the 
mean values. The independence of fLBGa on the chosen best fit 
model clearly shows that the effects of fesc, ( Xhi ) and f ct /f c 
are degenerate on fLBGa- 


large rates of star formation, leading to Muv values that 
scale with M,. Secondly, we find that both fLBG and fLBGa 
increase with an increase in the UV magnitude (or M*), as 
explained in Sec. |4.1| above: fLBG increases from 65% to 
100% as M* increases from 10 8 * — 10 10 ' 5 Mq. As a result of 
the stricter luminosity and EW criterion imposed to identify 
galaxies as LAEs, fLBGa < fLBG and declines more rapidly 
than fLBG towards fainter UV luminosities, decreasing from 
100% to 10-30% as M* decreases from 10 10 ' 5 — 1O S M 0 . Fur¬ 
ther, the decrease in L™ 1 with increasing / esc results in a 
(linear) decrease in /lab as shown from panels (a)-(c) of the 
same figure: while /lab is comparable for f eac = 0.05, 0.25, 
it decreases by about 0.15 for f esc = 0.5, where half of the 
ionizing photons do not contribute to the Lya luminosity 
thereby reducing the fraction of time it shows Lya emission. 

We then calculate /lbg and /lbgh including the effects 
of dust and IGM attenuation for each of the best-fit parame¬ 
ter combinations that match both the observed Lya LF and 


ACF as shown in Sec. [3] We use the f c value for each galaxy 
according to its final dust mass at z ~ 6.6, the Lya trans¬ 
mission T a of each galaxy was obtained from the ionization 
field and the ratio of the escape fractions of Lya and UV 
continuum photons (fa/fc) was set according to Table[l] As 
can be seen from Fig. [5] (since fLBGa is nearly identical for 
all best fit cases, we show only one best fit case) the addi¬ 
tional attenuation by dust in the ISM and neutral hydrogen 
in the IGM leads to a rise of the mean stellar mass in each 
Muv bin, as well as to lower values for /lbGq (compared to 
the intrinsic case) but not for /lbg- Nevertheless, we find 
the same trends as when considering intrinsic luminosities: 
/lbg always exceeds fLBGa , and /lbg and /lbgq decrease 
towards fainter UV luminosities with fLBGa declining more 
rapidly. However the relative difference between fLBG and 
fLBGa is larger and the decline of fLBGa is more rapid: while 
fLBG decreases from 100% to 65% as Muv increases from 
-23.5 to -18.5, fLBGa drops from 100% to essentially 0 for 
the same magnitude range. The stronger decline in fLBGa 
is not only due to the additional dust- and IGM attenua¬ 
tion of the Lya luminosities, but also due to the additional 
selection criterion in Lya equivalent width which becomes 
more important towards UV fainter galaxies. 

We also find that the decrease in f esc is compensated by 
the increase in fa/fc-, leading to very similar fLBG,a ratios 
for all the models; note that the IGM is almost ionized in 
most cases, leading to similar T a values. This clearly shows 
that (xhi), fesc and f a /fc compensate each other (as shown 
Hutter et ah||2014|, as a result of which fLBGa only de¬ 


pends on the combination of the parameters but not on their 
individual values. 

Thus in our model, the most luminous (massive) LBGs 
most often show Lya emission, irrespective of whether in¬ 
trinsic or observed luminosities are considered. 


5 DISCUSSION & CONCLUSIONS 


We couple a cosmological hydrodynamical simulation 
(GADGET-2) with a dust model and a radiative transfer 
code (pCRASH) to model high-z LAEs. Starting from a neu¬ 
tral IGM, we run pCRASH until the IGM is completely ion¬ 
ized, for fesc values ranging from 0.05 to 0.95. In [Hutter et al.| 
12014), we showed that comparing model results to Lya LF 


observations simultaneously constrains the escape fraction of 
ionizing photons fesc, the mean amount of neutral hydrogen 
( xhi) and the ratio of the escape fractions of Lya photons 
and UV continuum photons f a /fc to (xhi }— 1CU 4 — 0.5, 
fesc = 0.05 — 0.5 and f a /fc = 0.6 — 1.8. In this paper, 
we calculate the ACFs for these different combinations and 
find that comparing these to observations significantly nar¬ 
rows the allowed 3D parameter space (within a 3<r error) to 
(Xhi)— 0.01 - 10- 4 , fesc = 0.05 - 0.5 and f a /f c = 0.6 - 1.2. 
While the effects of these three parameters are degenerate 
on the Lya LFs, the ACF is most sensitive to large-scale 
ionization topologies and reionization leaves clearly distin¬ 
guishable ACF imprints (boosting up the strength of the 
ACF) that can not be compensated by varying / esc or f a /fc- 
Further, the ACF allows us to constrain ( Xhi )?S 0.01, inde- 
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Figure 6. Fraction of LBGs detected also as LAE (Xlab) as 
a function of UV luminosity. The fractions are shown for fiducial 
best fit models of fesc = 0.05: ( Xhi)= 0.01, fa/fc = 0.68 (red 
dashed) and ( Xhi)= 10 4 , fa/fc = 0.60 (blue solid). Error bars 
are derived from the Poissonian errors of the LAE and LBG num¬ 
bers. The upper panel (a) shows the fraction Xlae for a LAE 
selection criterion of L a > 10 42 erg s _1 and EW > 20A, and 
the lower panel (b) adopts an selection criterion of EW > 55A, 
which is in agreement with the observations by |Stark et a.l.|( 201 Oj). 
Black points represent the observational constraints of|Stark et al.| 
poTo) at 4.5 < 2 < 6.0. 


pendent of the other two parameters, and we also constrain 
/esc ^ 0.5 and f a /f c ^ 1.2. 

We then analyse the average time evolution of UV 
and Ly a luminosities of simulated galaxies at z ~ 6.6 in 
three bins of M* = 10 8 ’ 9 ’ 10 Mq, finding the following: as 
soon as a galaxy exceeds a critical stellar mass of M* ~ 
10 8 ' 5 (10' ,5 )Mq, its intrinsic Lya (UV) luminosity is large 
enough for it to be identified as a LAE (LBG). Including 
the effects of dust and IGM attenuation naturally results in 
an increase in this critical mass to 10 9,5 and 10 8 ' 5 Mq for 
LAEs and LBGs respectively. 

Considering the fraction of time during the last 
100 Myrs (prior to 2 = 6.6) a galaxy spends as a LBG with 
Lya emission ( fLBGa ) or as a LBG (/lbg), we find that the 
former is always smaller due to the more stringent luminos- 
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ity and EW constraints imposed on identifying galaxies as 
LAEs. We find that both the intrinsic and dust-attenuated 
fraction fLBGa and /l_bg rise with increasing UV luminos¬ 
ity (and hence A/*): intrinsically, /lbg (fLBGa) increases 
from 65% to 100% (10-30% to 100%) as M , increases from 
10 8 — 10 1o,5 Mq. As expected, including the effects of dust 
and IGM transmission reduces the values for fLBGa such 
that fLBGa decreases from 100% to essentially 0 as M* de¬ 
creases from 10 10 ' 5 — 10 8 Mq. Finally, we find that the frac¬ 
tion fLBGa of all our models that reproduce the observed 
Lya LF and LAE ACFs are independent of the chosen set of 
parameters: a larger f a /fc compensates a decrease in T a , or 
an increase in f eS c- As a result, fLBGa only depends on the 
combination of these 3 parameters but not on their individ¬ 
ual values. Thus, it is most often the most luminous LBGs 
that are visible in the Lya . 


Finally, we summarise the major caveats involved in 
this study. Firstly, given the cosmological volumes probed 
by the simulation, we are unable to resolve Lyman Limit 
systems (LLS) which could lead to a further decrease in the 
transmission T a along lines of sight (LOS) intercepted by 
such systems (Bolton & Hachnclt 2013). However, whether 
LLS are preferentially located in clustered regions, leading to 
an increasing suppression of T a for massive galaxies remains 
an open question. 


Secondly, as a natural consequence of simulating cosmo¬ 
logical volumes we are unable to resolve the ISM of individ¬ 
ual galaxies, for which reason we assume a Gaussian profile 
(with a width set by the rotation velocity of the galaxy) for 
the Lya line that emerges out of any galaxy which is proba¬ 
bly an unrealistic scenario (see e.g. Verhamme et al.] p 068 |. 
We note that our constraint on the ionization state of the 
IGM is model dependent, since the IGM Lya transmission is 
sensitive to the the assumed line profile (Jensen et al.|2013|. 


Thirdly, we assume dust attenuation and IGM trans¬ 
mission to be equal to the values at z ~ 6.6 in order to 
calculate both fLBGa and /ubg- While the dust mass (and 
hence attenuation) would be expected to be lower at earlier 
times, tracing this buildup would require tracking the dust 
growth in the progenitors of our simulated galaxies. This is 
beyond the scope of our present paper and we defer to this 
analysis to a work that is in preparation. Fourthly, fLBGa 
would be expected to decrease with increasing z as a re¬ 
sult of an increase in (xhi) (leading to a decrease in T a ). 
However, properly accounting for the latter effect requires 
modelling the entire history of reionization. 


While we have explored the full range of possible values 
for f eac , its mass and 2 -dependence remain poorly under¬ 
stood, which is also one of the main caveats involved in mod¬ 
elling the time-evolution of reionization. An increase in f e3C 
with decreasing mass (Ferrara & Loeb 2013 e.g.) would sup¬ 
press the visibility of low-mass objects, and strongly impact 
the reionization fields we generate, emphasising the strong 


clustering of high-mass halos (Kaiser 1984 Bardeen et al. 


1986 Mo & White 1996), whilst depressing the LF at the 


faint end. Bringing these values into agreement with ob¬ 
servations would then require an f a /fc ratio that decreases 
with increasing mass. With its observations of the ionization 
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topology, instruments such as LoFAR will be invaluable in 
answering some of these outstanding questions. 

Finally, we calculate the fraction of LBGs that would 
be identified as LAEs, Xlae ■ Imposing a LAE selection cri¬ 
terion of L a ^ 10 42 erg s _1 and EW ^ 20A we find that it is 



Xlae does not show the behaviour observed by [Stark et al.| 
(2010), increasing instead of decreasing with UV magnitude. 
This mis-match is probably due to a combination of the fol¬ 
lowing effects: firstly, the Lya IGM transmission ( T a ) is sub¬ 
ject to a large variance along different lines of sight (LOS) 
(see Hutter et al. [20141 due to the patchy nature of reion¬ 
ization; while we use T a values averaged over 48 LOS, us¬ 
ing values along a specific LOS would lead to an over- or 
underestimate of T a . Secondly, as mentioned before, it is 
possible that the inclusion of LLS could decrease T a of mas¬ 
sive galaxies, leading to lower Lya EW values respectively. 
Thirdly, the results by Stark et al. (2010) are based on a 
post-reionization sample of galaxies (4.5 < z < 6.0), while 
our simulation samples the end of reionization era. The evo¬ 
lution of the UV (and Lya ) LFs of LAEs and LBGs sug¬ 
gest that the emitted observable Lya radiation varies with 
cosmic epoch, depending predominantly on the evolution of 
dust and gas at « < 6 and on the IGM neutral hydrogen 
fraction at 2 > 6. Thus, the probe of LBGs with/without 
Lya emission in Stark et al. (2010) (dusty, ionized IGM) 
may differ to the LBG and LAE population at the end of 
reionization (less dusty, partly neutral IGM). Indeed, a re¬ 
markable high fraction of strong LAEs among 6.0 < 2 < 6.5 
luminous LBGs (—21.75 < Muv < —20.25) (|Curtis-Lake 


|e~ al.| |2012| indicates that most LBGs are also identified 
as LAEs, which is in agreement with our findings. How¬ 
ever, their sample suffers from low statistics, which will be 
overcome by upcoming Subaru/HST/UltraVISTA data. Fi¬ 
nally, as an observational caveat, Verhamme et al. (2012) 
have shown that detailed Lya radiative transfer calculations 
of simulated galaxies suggest stronger inclination effects for 
Lya photons than for UV continuum photons, introducing 
biases in the selection function of narrow-band LAE surveys 
that could lead to a significant fraction of LBGs galaxies be¬ 
ing missed as LAEs. We aim at investigating these effects in 
detail in future works and shedding light on the tantalising 
connection between LBGs visible/invisible as LAEs. 
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APPENDIX A: THE FRACTION OF LIFETIME 
SPENT AS LBG WITH AND WITHOUT LYa 
EMISSION 


We also show results for the limiting case where the age is 
the time since the onset of star formation, showing how the 
fractions of lifetime change when the time span is increased. 
We take the total age of a galaxy as the time between the 
formation of the first star and z = 6.6, and calculate the 
corresponding fractions of lifetime that galaxies in different 
Muv bins spend as a LBG and as a LBG with Lya emission 
If lbg and fLBGot , respectively). As described in Section 
|4.2| we compute the fractions of lifetime for our best fit 
models (see Table [lj and show the fractions for the f esc = 
0-5, (xhi)— 0.01 and f a /fc best fit case in Fig. Al We 
note that the fractions fLBGa of all other best fit cases are 
nearly identical to the shown one. Compared to the fractions 
derived for a time span of 100 Myrs, the fractions of the total 
lifetime decrease for all UV luminosities. Galaxies that have 
been identified as LBGs or LBGs with Lya emission within 
100 Myrs prior to z = 6.6 have been fainter at earlier times 
where they have not met the selection criteria. The decrease 
in fLBG and fLBGa illustrates the growth of stellar mass 
in galaxies with time and the existence of a critical stellar- 
mass above which a galaxy produces enough luminosity to 
be identified as a LBG or as a LBG with Lya emission. 
However, we like to note that the point in time when the 
first star forms depends on the resolution of the simulation, 
making our fractions of lifetime resolution-dependent but 
showing their trends as the time span is increased. 
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Figure Al. Fraction of lifetime that galaxies spend as LBGs 
(fLBG 5 circles) or as LBGs and LAEs (fLBGa ? squares) as a 
function of the UV luminosity for our best fit models. The mean 
stellar mass in each Muv bin is encoded in the shown colour scale. 
This panel shows the best fit case for f esc = 0.05, ( Xhi) = 10 —4 
and fa/fc = 0.60. We omit to show the other best fit cases, since 
they are nearly identical to the shown one. For each best fit model 
we assume individual f c values for each galaxy according to its 
final dust mass at z ~ 6.6; the Lya transmission T& of each galaxy 
was obtained from the respective ionization field and the ratio of 
the escape fractions of Lya and UV continuum photons was set 
according to Table [l] The fractions are computed as the mean of 
the galaxies within Muv bins k ranging from k — 0.25 to k + 0.25 
for k = —25 .. —18 in steps of 0.5. Error bars show the standard 
deviations of the mean values. 
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